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ABSTRACT 

Most stars form in a cluster environment. These stars are initially surrounded by discs from which 
potentially planetary systems form. Of all cluster environments starburst clusters are probably the 
most hostile for planetary systems in our Galaxy. The intense stellar radiation and extreme density 
favour rapid destruction of circumstellar discs via photoevaporation and stellar encounters. Evolving 
a virialized model of the Arches cluster in the Galactic tidal field we investigate the effect of stellar 
encounters on circumstellar discs in a prototypical starburst cluster. Despite its proximity to the 
deep gravitational potential of the Galactic centre only a moderate fraction of members escapes to 
form an extended pair of tidal tails. Our simulations show that encounters destroy one third of the 
circumstellar discs in the cluster core within the first 2.5 Myr of evolution, preferentially affecting the 
least and most massive stars. A small fraction of these events causes rapid ejection and the formation 
of a weaker second pair of tidal tails that is overpopulatcd by disc-poor stars. Two predictions arise 
from our study: (i) If not destroyed by photoevaporation protoplanetary discs of massive late B- and 
early 0-type stars represent the most likely hosts of planet formation in starburst clusters, (ii) Multi- 
epoch K- and L-band photometry of the Arches cluster would provide the kinematically selected 
membership sample required to detect the additional pair of disc-poor tidal tails. 
Keywords: methods: numerical - stars: kinematics and dynamics, pre-main sequence 



1. INTRODUGTION 

Observations in the past decade have shown that most 
young stars do not for m in isolation but as part of a clus- 
ter e nvironment (e.g. |Lada fc Lada 2003 Porras et al. 

103 Evans et al. 2069p . I'he accretion discs ot these 
stars are thus exposed to environmental effects that could 
affect their evolution (iR ichling & Yorkc 1998; Hollcnbach] 
et al.||2000 ^, 'Scally & Clarke 2001 ; Ada ms et al.l | 20d^ 
Olczak et al. 2006; Pfalzner 2006; Pfalzner et al.ll2006 
F!alzncr & Olczak 2007; Clarke 2007; Olczak ct al. 2010) 
A^s these discs are the prerequisites for the formation of 
planetary systems this process might be influenced by the 
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cluster environment. However, so far theoretical investi- 
gations about the effect of irradiation by massive stars 
or strong gravitational interactions have concentrated on 
low- and intermediate mass star clusters like IC348 or the 
Orion Nebula Cluster (ONC). Much more massive sys- 
tems like NGC 3603, the Arches cluster or Wester lund 1 
- known as starburst clusters - that are expected to trig- 
ger the strongest effects have not been treated. Con- 
sidering the huge amount of massive stars and extreme 
densities, i.e. more than 100 0-stars an d a core d ensity 
>10^Mopc~^ in the Arches cluster (^Figer et al.|[T999|), 



an extrapolation of the previous investigations towards 
starburst clusters would suggest that the lifetime of cir- 
cumstellar discs could be shortened dramatically and so 
planet formation in such an environment hindered signif- 
icantly. 

However, recent observations of the Arches clus ter in 
the near-infrared JKL' bands by Stolte et al. ( 2010 1 have 
revealed emission from circumstellar matter around at 
least a few percent of stars in the mass range 2 — 20 M0 . 
This detection of discs in the Arches B-star population 



2 



Olczak et al. 



was surprising for two reasons. First, it is expected 
that the stars' UV radiation causes depletion of their 
own inner disc. Considering that the characteristic UV 
evaporation timescale of a pri mordial disc around Her - 
big Be stars is less than 1 Myr (|Alonso-Albi et"aL]|2009[ ), 
a disc lifetime of 2.5 Myr for B-type stars implies that 
the self-photoerosion of discs is probably less efficient 
than various mo dels suggest (e.g. Hollenbach et al.|2000 



Alexander 2008). Second, extrapolating from the above 
mentioned numerical studies of less massive clusters, in 
a starburst cluster environment disc destruction is ex- 
pected to be much increased by external processes. The 
extreme UV radiation from numerous 0-stars and grav- 
itational interactions induced by the high stellar density 
could boost the removal of external disc material. 

In the present investigation we focus on the mecha- 
nism of encounter-induced disc-destruction to constrain 
its contribution to the overall disc life time. We use sim- 
ulations of the Arche s cluster preformed by the starlab 
simulation package ( McMillan 1996 Portegies Zwart^ 
et al. 2001 Hut 2003p and carry out the analysis of 
the disc-mass loss analogous to our previous publications 
referenced above. Hence this study of a massive clus- 
ter complements previous investigations of encounter- 
induced disc evolution in low- and intermediate-mass 



clusters (jScally k C larke 2001 [Adams et a l. 2004; Olczak 
et al.||200^! 'Pfalzner 2006; Pt'alz ner et aTl2006t jPfalzner 



czak 2007; Clarke 2007; Olczak et al.||2010| 
In contrast to the low- and intermediate-mass clus 
ters in the solar neighbourhood the location of the mas- 
sive Arches cluster near the Galactic centre exposes it 
to strong tidal fields. Therefore another difference to 
previous simulations is that our model includes the con- 
tribution of the Galactic tidal field. Stars escaping from 
a star cluster in the gravitational field of a galaxy form 
extended tidal t ails known from various observations in 
the Milky Way (iGril lmair et al.lll995MKharc henko et al 
19971 iLeon et al. 2DD0; Belokurov et al. 2006) and numer- 



ical models (ICo rnFes et al. 1999; Porteg ies Zwart et al . 
2002||Capuzzo Dolcetta et al. 2005; FcllhauereFaT 



Except for high-speed escapers created in Tcw-body en- 
counters, stars escape as a result of two-body encounters 
passing at slow speed close to the saddle point s of the 
effective po tential, known as Lagrange points (Kiipper 
eraL][2008l ). 



With our study we aim to investigate the evolution 
of the encounter induced-disc mass loss in a tidally dis- 
torted starburst cluster and its tidal arms. Of particular 
relevance is the work of Just et al. ( 2009 ) who use the 



epicycle theory for a quantitative analysis of the tidal tail 
structure of star clusters moving on a circular orbit in the 
Galactic disc. They find that the radial offset of a star's 
epicyclic motion relative to the orbit of the star cluster, 
Ai?o, is first order in the angular momentum, AL, while 
its radial amplitude rm depends on the energy excess, 
Ai?, which is of second order in AL: 



Ai?o oc AL , 

Tni OC VAi?. 



(1) 



Just et al. 



( 2009[ ) conclude that the circumstances are 
more complicated near the Galactic Centre but to lowest 
order the theory is still applicable. 
Throughout this work we assume that initially all stars 



are surrounded by protoplanetary discs. This is justi- 
fied by observations that reveal disc fract ions of nearly 
100% in very young star clusters (e.g. Haisch et al. 
20001 |Lada et al.|[2000l [Haisch et aL][200l| 1Hillenbrana' 
2005 1 . The typical disc diameter is a few hundred AU 



for low- and intermediate-mass sta rs (McCaughrcan fc] 
0'dell[[1996| [Andrews fc Wimams[ |2007| ) , though discs 
ot more than several thousa nd A U surrou nding massive 
stars have been observed (see Zhang|20^ and references 
therein). However, it remains unclear whether a clear 
correlati on of disc extension an d stellar mass does exist 
(see e.g. [Vic ente fc Alves[|2005| . 

In Section |5] we outline the observationally determined 
basic properties of the Arches cluster that serves as a 
reference for our cluster models. The computational 
method and the properties of the numerical models are 
described in Section [3] Afterwards we present results 
from our numerical simulations in Section [4| The con- 
clusion and discussion mark the last section of this paper. 

2. PROPERTIES OF THE ARCHES CLUSTER 

The Arches cluster is one of the most massive and dens- 
est youn g clusters i n the Galaxy and harbors about 125 
O stars (Figer et al. 1999). Most of its stars have formed 



more or less simultaneously abou t 2.5 Myr ago (Najarro 
et al.[[2004] iMartins et al.[|2008[ ). The cluster is very 
comp act witn a core radius » 0.15 pcjEspinoza et al 
2009|), a half-mas s radius i?hm ~ 0.4 pc (jFig er ct al. 2002" 
Stolte et al.|2005[) and a tidal radi us Rt ^, 1 pc ( Kim et at 

(2006* 



Bin- 



Based on these size scales Espinoza et al 
quote a concentration parameter c = log [Kt/Kc) 

This is equivalent to a King parameter Wq ~ 4 (e.g 

ney fc Tremaine[1987 ). In their detailed analysis Hartst 
et al. " pblo| ) find that stars with m > 10 Mq are more 
concentrated due to mass segregation and best fit by a 
King profile with Wq = 7. 

The Arches cluster is located at a proje cted distance 
to the Galactic Center of only about 30 pc ( Nagata et al. 
19951. Taking into accou nt its three-dimens ional space 



motion of 232 ± 30km/s (Stolte et al. 2008), numerical 
simulations suggest an eccen tric orbit and a true present- 
day distance below 200 pc (Stolte et al. 2008; see also 
Portegies Zwart et al. 2002^ ! 'I'he internal motion of 



the cluster is characterized b y a one-dimensional v elo city 
dispersion of 5.4±0.4km/s ( |clarkson et"aI][20TT| ). Note 
that these authors have refined the bulk proper motion 
of the Arches cluster that implies a slightly smaller space 
motion of 196 ± 17km/s. 

It is still under debate whether there is an overabun- 
dance of massive stars, i.e. an unusually flat ("top- 
heavy") mass function, or if the mass function is trun- 
cated at_jnasseswen_above the hydrogen burning limit 
(e.g. Stolt e et al. 2005 ). However, in their most re- 



cent work Espinoza et al.| ( [2009) do not find any sig- 
nificant difference to the less massive star-forming re- 
gions in the solar neighbourhood. They derive a mass of 
about 3 X 10^ McT) within the estimated half-mass radius 



(assuming the IMF of [Kroupa 2001) and a core density 
of 2 X lO^M0pc~^. However, these estimates rely on 
counts of stellar masses above 10 Mq only. The same 
authors find evidence for mass segregation in the clus- 



ter (e.g. Stolte et al. 2005 Espinoza et al. 20091. De 



spite its young age tfiis feature could be driv en dynami 



cally ( ,Spitzer„1969i |Farouki fc Salpeter[[1982[ ) due to the 
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Arches cluster' s short half-mass relaxation time of only 



2006 Portegies Zwart et al. 



, though t here is no consen sus 



effect (e.g. |Dib et al.||2007 1 



about lOMyr ( Kim et al 
2007l|Espinoza et al.|20()^ 

on the signi hcance or the e 

Recently, Stolte et al. ( 2010| ) determmed proper mo- 
tions of hundreds of stars and simultaneously obtained 
deep L'-band photometry of OB(A)-type stars for a pre- 
cise and reliable estimate of the disc fraction. They de- 
tect 21 sources in the mass range ~ 2-20 M0 with emis- 
sion in excess of the expected photospheric flux, indi- 
cating the presence of circumstellar dust, and find evi- 
dence for three optically thick discs from CO band head 
emission. This finding implies a total disc fraction of 
6 % ± 2 % within the observed mass range. The disc frac- 
tion increases with cluster radius from 2.7% ± 1.8% in 
the core (r < 0.16 pc) to 5.4% ± 2.6% at intermediate 
radn (0.16 pc < r < 0.3 pc) and 9.7% ± 3.7% outside 
(r > 0.3pc). 

3. COMPUTATIONAL METHOD 
3.1. Dynamical cluster model 

The dynamical mod el of the Arches clu ster used here is 
based on the work of Harfst et al. (2010) who compared 
in detail the results of numerical simulations with obser- 
vational data. They were able to constrain the initial 
conditions and found a best-fitting model as a function 
of five initial parameters: (i) mass function slope a, (ii) 
lower mass limit miow, (hi) number of massive stars A^ms 
(with mass m > 20 Mq), (iv) virial radius i?vir: a-nd (v) 
King parameter Wq. We summarize the main properties 
of this model. 

We assume the cluster is initially (at t — OMyr) gas- 
free and in virial equilibrium (i.e. the ratio of kinetic 
to potential energy is Qvir = 0.5) and t hat its ma ss is 
distributed according to a King profile (King 1966). A 
single-aged stellar population is used and no initial mass 
segregation is taken into account. Primordial binaries are 
not included but we discuss their potential dynamical ef- 
fects in Section [5] A grid of numerical models with vary- 
ing initial parameters a, miow, AmS: -Rvir, and Wq was 
evolved until the present-day age of the Arches cluster, 
^Ac = 2.5 Myr, to find the best match with observational 
data. 



Acc ording to the best-fitting models of Harfst et al. 
(2010) the in itial slop e a abo ve 0.5 M0 is consistent with 
the IMF of |Kroupa| ( |2001[ ). The models also favour 
a narrow range tor the initial number of massive stars 
Ams ~ 150 — 200. However, the initial lower mass 
limit miow is not well constrained by their models and 
thus only allows to estimate a lower limit for the initial 
total mass Af > 4 • 10*^ M©. The initial size of the cluster 
with i?vir « 0.7 — 0.8 pc is again well constrained. The 
corresponding core radius is about 0.4 pc initially, and 
shrinks, due the dynamical evolution of the cluster, down 
to the observed 0.2 pc. The best-fitting models have 
modest initial concentrations with Wq = 3 — 5. In sum- 
mary, the favoured fiducial model of the Arches cluster, 
which we will use here, is parametrized by Ams — 150, 
i?vir = 0.7, and Wo = 3. The total number of particles 
is A = 74153. 

Stellar evolution and the orbit of the Arches cluster in 
the Galactic cente r potentia l have been neglected in the 
simulations of .Harfst et al. (2010). Here we extend the 



fiducial model to lower stellar masses (miow — 0.1 M©) 
and include these two effects to generate two models: 
an "isolated" model I with stellar evolution only and an 
"orbital" model O with both stellar evolution and a real- 
istic orbital motion of the cluster in the Galactic center 
potential. 

The simulations were carried out with the direct N- 



Ian 


1996 ] 


mc 


udes m 



Portegies Zwart et al.||2001[ |Hut||20Ci3l ) that 
lodules for stellar evolution and an external 



stellar evolution and an external 
gravitational potential. The stellar evolution module ac- 
counts for mass loss by stellar winds and binary evo- 



lution dEggleton et"ar] |1989UPortegies Zwart fc Verbunt 
1996||Tout etal.|l997l|La 



tegL 



anger 199^r^ 'he external poten 
tial that was used in our simulations is that of a power 
law mass function with Me-ai (r) = 4.25 • I0^{r[pc])'^-'^ Ma 



( Portegies Zwart et al.| |2002) which is base d on observa- 
tions o t the A -band luminosity function by |Mezger et al.| 
( 1999 1 . Computations were acce lerated using CPUs w ith 
the help of the Sapporo library ( Gaburov et al.|[2009 ). 

The cluster orbit in the potential is given by the six 
phase-space coordinates at a given time. Five present- 
day (^Ac = 2.5 Myr) coordinates of the Arches cluster 
are well kno wn from observati ons: its line-of-sight veloc- 
ity (95km/s: Figer et al.|2002 ), prop er motion (19 0km/s 
anti-parallel to the Galactic plane: Stolte et al.| [2008), 
and projected posit ion (30 pc from the Galactic center: 
Nagata et al. 1995). We define a coordinate system in 
which the a:-axis is along the Galactic plane, the y-axis 
along the line-of-sight, and the 2-axis towards the Galac- 
tic north pole, to obtain the following present-day posi- 
tion and velocity vectors of the cluster: 



^•cluster = (-24, dios, 10) pC , 

i^ciustor = (-190, 95, 0) km/s , 



(2) 



where we adopt d\os — —100 pc as a probable value for the 
line-of- sight distance of th e Arches cluster to the Galactic 
centre ( Stolte et al.|2008 ). Throughout the paper we use 
Rna = okpc tor the Sun s distance to th e Galactic center 
( |Ghez et al.|[2008| [Gillessen et al][2009l ). 

'I'he present-day phase space coordinates were used to 
numerically integerate the orbit backwards in time to 
find the initial position and velocity of the Arches cluster. 
From there, the full cluster was then integrated with kira 
for 6 Myr including the effects from stellar evolution. Its 
orbital motion is shown in Fig. [T] 

We point out that model I is the best-fitting initial 
model of the present-day Arches cluster from our set of 
isolated models. However, computational limitations did 
not allow to carry out a second parameter study of orbital 
models. Hence model O is just a modification of model I 
that is evolved on a realistic orbit around the Galactic 
Centre in the Galactic tidal field. As we will show later 



in Section 4.1 this modification has a strong effect on 
the internal dynamics of model O such that it becomes 
a less good model of the internal dynamics of the Arches 
cluster compared to model I. 

3.2. Star-disc encounters 

The key aspect of this work is the determination of 
the encounter-induced disc-mass loss in a starburst clus- 
ter environment. For this purpose two different types of 
numerical simulations have been combined. First, as part 
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Figure 1. Three planar projections of the orbital motion of the cluster's centre-of-mass for model O. The origin marks the Galactic centre. 
The grey scale and the line thickness decode the time evolution over 6Myr; the black/grey filled circles represent the initial/final position, 
respectively. The position at the present-day age of the Arches cluster (Sac = 2.5 Myr) is indicated by the dark grey filled rectangle. 

of previous work a parametrized fit of the relative disc- 
mass loss has been derived from a tree-code based pa- 



rameter study of isolated star-disc e ncounters ( Pfalzner 
et al.||2005[ |2006[ [Olczak et al |J2010[ ). Using this tit for- 
mula, m a second step the individual relative disc-mass 
loss of star cluster members was determined from direct 
./V-body simulations of star cluster dynamics, in which 
the encounter history of all stars was tracked. We out- 
line the procedure in more detail below. 

Stellar encounters in dense clusters can lead to signifi- 
cant transport o f mass and angular momentum in proto 
planetary discs (Olczak et al. 



parabolic derived by [Olczak et al.| f2010f : 

(5md(e) — Snid 

X exp[-0.12(e- 1)] 
X {0.83-0.015(£-l) 



Pfalzner fc Olczak.20070 nt; 



0.17exp[0.1(e- 1)]}. 

(4) 

Note that this fit is only a compromise: by using the 
median we underestimate the effect of strong perturba- 
tions (i.e. those induced by close appro aches of massive 
perturbers). However, as shown by, Olczak et al. (20101 



le present invest igation we 
have used Eg. (Ij from [Pfalzner et al. (20061 to deter- 
mine the encounter-induced relative disc-mass loss Sm^ 
of a star with mass Mi and disc radius caused by a 
perturber with mass M2 and pericentric distance Tpi 



2006[ [Pfalzner et al.||2006| this type of interactions is expected to occur rarely in 



M2 + 0.5Mi 



log 



2.8 1 



0.1 



X exp ■ 



Ml 

M2 + O.bMi 



3/2 



- 0.5 



(3) 



Note that Eq. (|3| is valid for low relative disc masses 
only, m<i < 10 Mi, implying negligible viscosity and 
self-gravity. We point out that in this limit the esti- 
mated relative disc-mass loss 6md is independent of the 
disc mass but depends on the disc radius. However, it 
represents an upper limit because the parameter study 
was restricted to co-planar, prograde, parabolic encoun- 
ters that are the most perturbing. We note that Eq. |3| 
has been derived for discs with a fixed surface density 
profile E (X 



Recently, Steinhausen et al. 



lensity 
( [2012| ) 



have shown that the shape of the disc-mass distribution 
has a significant impact on the quantity of the disc-mass 
and angular momentum losses in star-disc encounters. 
Maximum losses are generally obtained for initially flat 
distributed disc material. 

To account for the reduced disc-mass loss in hyperbolic 
encounters with eccentricity e > 1, we use the flt function 
of the median relative disc-mass loss, normalized to the 



stellar systems as dense as the Arches cluster. 

The pure stellar dynamical simulations of a cluster en- 
vironment have been carried out with the direct iV-body 
integrator kira as outlined in Sec. 3.1 The dynami- 



cal data of each particle (time t, mass Afi, position r, 
and velocity v) have been dumped every ten individual 
integration time steps to create a high-resolution tem- 
poral grid of particle phase space vectors. Finally, this 
data set was parsed by a dedicated encounter tracking 
software to generate for each particle a record of inter- 
action parameters with its strongest perturber (time t, 
particle mass Mi, perturber mass M2, separation r, and 
eccentricity e) at the beginning, the predicted pericentre 
passage, and the end of the interaction. The algorithmic 
basics of the procedure are described in detail in App.fX} 
In the remainder of the paper we will refer to normal- 
ized disc-mass rfid and normalized disc-mass loss Arhd- 
Normalizingthe disc-mass rhd € [0, 1] is justified when 
using Eqs. |3| and (Q that are independent of the abso- 
lute disc mass as outlined above. The normalized disc- 
mass loss Amd is defined as the absolute loss of the nor- 
malized disc-mass, 



A?77,(j = 777(1 ■ Sm^ 



(5) 



In our approach a simplified prescription assigns stars 
into one of two distinct groups: if the accumulated nor- 
malized disc-mass loss Arhd exceeds 0.9 (i.e. 90 % of 
the initial disc mass), stars are marked as "disc-less"; 
otherwise they are termed "star-disc systems" . Accord- 
ing to observations as outlined in Sec. [T] we fix the typ- 
ical disc diameter of a solar mass star to be 300 AU 
and consider two different model distributions of the ini- 
tial disc radii r^: i) scaling with stellar mass AIi as 
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Figure 2. Comparison of the dynamical evolution of models I 
(solid) and O (dotted). The black and grey lines mark the core 
density (left logarithmic scale) and the half-mass radius (right lin- 
ear scale), respectively. 



Td — 150 AUa/Mi/M0, which is equivalent to the as- 
sumption of an equal force at the disc boundary, and 
ii) fixed disc radii = 150 AU. Whenever results are 
presented, we will specify which of these two distribu- 
tions has been used. However, in the present study we 
focus on the latter model as no obvious correlation be- 
tween disc radius and stellar mass was found from ob- 

Note that 



2005D . 

; the distribution 



servations so far (see Vicente fc Alves 
we make no specific assumption abou 
of disc masses. As outlined before the normalized disc- 
mass loss Arhd is independent of the absolute disc mass 
in the limit of low relative disc masses, md < lO^^Mi, 
assumed here. 

4. RESULTS OF NUMERICAL SIMULATIONS 

For the two cluster models 1 and O a full mass distri- 
bution was gen erated according to the IMF published by 
Kroupa (2001 1. Using these two models we provide in 
the following a detailed view on the dynamical evolution 
of young stars and their circumstellar discs in the Arches 
cluster. 

4.1. Cluster dynamics 

The orbital model O is a refined version of the isolated 
model I that is evolved on a realistic orbit around the 
Galactic Centre in the Galactic tidal field. It serves to 
analyze the evolution of the tidal features of the Arches 
cluster. However, unlike model I, model O is not the best 
representation of the internal dynamics of the present- 
day Arches cluster because it is not in a dynamical equi- 
librium state initially. The imposed Galactic tidal field 
acts as an additional energy source such that model O 
expands initially and re-bounces slightly before setting 
into an equilibrium state. 

We demonstrate this effect in Fig. [2] Initially, the core 
density ricoro of model O (black dotted line) drops rapidly 
by a factor of three, then rises slightly and afterwards 
remains roughly a factor of two below the value of model I 
(solid black line) . The two lines are running parallel due 
to the logarithmic scale. An equivalent behaviour is seen 
for the half-mass radius (grey lines): the ratio of the 
half-mass radii remains at about 1.3 (that is a factor of 
two in volume) after the re-bounce. Note the divergence 



of the two lines due to the linear scale. Hence, after 
model O has reached dynamical equilibrium at ~0.3Myr 
the evolution of both models is dynamically equivalent 
over the remaining 5.7 Myr (or more than one hundred 
crossing times). 

In particular do both models share the same charac- 
teristic local maximum of the core density ncoro at the 
present-day cluster age ^ac ~ 2.5 Myr. This time tAC 
marks as well the minimum of the 1 % Lagrangian ra- 
dius and the formation of a massive binary. It is fol- 
lowed by a significant continuous expansion of the clus- 
ter as is evident from the evolution of the half-mass 
radius. Thes e features are characteristic signatures of 
core collapse ( Lynden-Bell fc Wood||1968 Farouki et al. 
1983 1 . We thus identity the expansion of the cluster after 
~2.5 Myr of dynamical evolution as its post-core-collapse 
phase. 

4.2. Encounter dynamics 

The main purpose of the present work is a realistic 
study of the encounter-induced disc-mass loss of stars 
in the Arches cluster by taking into account the exter- 
nal effect of the Galactic tidal field. However, we have 
demonstrated in the previous section that the dynamical 
model with tidal field, model O, is too extended com- 
pared to model I, our best-fitting isolated model of the 
observed properties of the Arches cluster. Consequently, 
the encounter rate of star-disc systems in the underdense 
model O is expected to be too low. Ignoring for a sim- 
plified dimensional argument the effect of gravitational 
focusing, the frequency of encounters with distance r de- 
creases with increasing radius i? of a self-gravitating sys- 
tem, 

(6) 



fr = TicrS oc R-^R-^f^r^ = R-^'^r^ 



where n, cr, and S denote its density, velocity dispersion, 
and cross section with radius r. The nearly constant size 
ratio of models O and I over time, TZ = R'~^ /R^ ~ 1.3, 
implies a roughly time-invariant ratio of encounter fre- 
quencies F = f^/ff (X 7^-^/2 « 0.4. However, Eq. ^ 
shows that we can compensate for this difference by re- 
ducing the encounter distances r in the calculation of the 
disc-mass loss for model O by a factor J^^/^ 0.6. In 
fact, technically we have increased the standard disc ra- 



3.2 



dius for model O by to r^ = 250 AU (see Sec. 

for reference) and found a very good agreement between 
the encounter histories of models O and I. 

As will be shown later the only significant difference in 
the effect of encounters between models O and I that can 
not be corrected for is gravitational focusing. Its strong 
dependence on stellar mass implies that the encounter- 
induced disc-mass loss of massive stars is reduced in 
model O. However, the general results for the entire 
cluster population are not much affected by the rela- 
tively rare interactions of high-mass stars as discussed 
in Sec. [121 

The overall effect of the encounter-induced disc-mass 
loss in the Arches cluster is shown in Fig. [3] Here the 
temporal evolution of the cluster disc fraction (CDF) is 
plotted. Note that stars with a normalized disc-mass loss 
of more than 0.9 are defined as disc-less. The upper plot 
has been calculated using Eq. ^ , restricted to parabolic 
encounters, while the bottom plot represents a more re- 
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Figure 3. Time evolution of tiie encounter-driven cluster disc 
fraction (CDF) in two different volumes of model O accounting for 
a) parabolic encounters via Eq. ijsjl and b) eccentric encounters 
via Eq. (j4]l. Grey lines mark the core and black lines two times the 
half-mass radius at a given time. Const a.nt i nitial disc radii have 
been assumed for this calculation (see Sec. |3.2[ l. The vertical dotted 
line marks the estimated age of the Arches cluster. The time scale 
is shown in both physical (lower scale) and Af-body units (upper 
scale). 

alistic scenario considering the eccentricity of encounters 
via Eq. Here constant initial disc radii were assumed. 

From Vig. [3]a one would expect that more than 80 % 
of all discs are destroyed after 2.5 Myr, and after 5Myr 
basically no discs would have survived the repeating ero- 
sion by parabolic gravitational interactions. This sce- 
nario suggests that stellar encounters are potentially a 
very efficient disc destruction mechanism in starburst 
clusters that could suppress planet formation entirely. 
However, unlike for intermediate-mass clusters like the 
ONC the assumption of parabolic encounters is a strong 
simplificati on for dense stellar systems like the Arches 



cluster (see Olczak et ah] 2010) and thus largely overes- 
timates the disruptive ettect on discs in this cluster. In 
hyperbolic encounters the interaction time between the 
encounter partners is shorter and thus reduces t he trans- 
port of energy and angular momentum (see e.g. Ostriker 
19941 ). 

Wnen we take the eccentricity of the stellar encounters 
into account (Fig. Isb) the effect on the disc-mass loss is 
much weaker, resulting in the destruction of 10 % of the 
discs in the entire cluster (here defined as two times the 
half-mass radius) and roughly 30 % in the core up to its 




eccentricity 

Figure 4. Number of encounters as a function of eccentricity 
for model O. Constant initia l disc radii have been assumed for 
this calculation (see Sec. |3.2[ ). The different shaded regions mark 
all 516243 encounters (light grey), 31044 encounters causing more 
than 50 % (grey) and 2412 encounters with more than 90 % disc- 
mass loss (dark grey). 



present age. However, not only in the core but in the 
entire cluster the encounter-induced destruction of cir- 
cumstellar discs does continue over the entire simulation 
time and leads to a fractional disc destruction of 25% 
in the entire cluster and 50 % in the core after 6 Myr of 
dynamical evolution. 

The characteristic onset of expansion at ^ 2.5 Myr as 
identified in Fig. [2] manifests itself in the kink of the 
core disc fraction in the bottom panel of Fig. |3] (grey 
line). The strong decrease in density reduces the stellar 
encounter rate and thus the rate of disc destruction. 

We find that encounters of young star-disc systems in 
the Arches cluster can only destroy a minor fraction of all 
protoplanetary discs within a few Myr. The encounter- 
induced disc-mass loss is drastically reduced by a large 
fraction of highly eccentric fly-bys. The reason is that hy- 
perbolic encounters typically result from chance encoun- 
ters of low- and intermediate-mass stars while parabolic 
encounters are typically the outcome of significant gravi 
tational focusin g of low-mass stars by massive stars ( 01 



czak et al. 2010 ). The corresponding encounter rates jhyp 
and /par iiave very different scaling behaviour with the 
cluster density p: /hyp oc p'^/^ and /par oc p^/^. While 
the intermediate-mass ONC is dynamically balance d by 
hyperbolic and parabolic encounters (see Fig. 5 of 01- 
czak et aL||2010 ), in the roughly 20 times denser Arches 
cluster the traction of parabolic encounters is less than 
10%. Hence the relative importance of gravitational fo- 
cusing decreases very rapidly with increasing cluster den- 
sity such that hyperbolic encounters become by far the 
dominant type of stellar interactions in a starburst clus- 
ter. This is demonstrated for our model O of the Arches 
cluster in Fig. |4j Here we plot the number of encounters 
as a function of eccentricity for three different thresholds 
on the normalized disc-mass loss Amd (marked by dif- 
ferent grey scales) . About 73 % of all encounters (light 
grey) have large eccentricities e > 10. 

The following plots are thus based on calculations that 
take the eccentricity of the stellar encounters explicitly 
into account. In Fig. [5] we show the average normal- 
ized disc-mass loss as afunction of stellar mass for two 
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b)i 



mass [Msu„] 

Figure 5. Average normalized disc-mass loss as a function of 
stellar mass after 2.5 Myr of simulation time for a) model I and 
b) model O. The black and grey bars mark the sce nario s using 
constant and scaled initial disc radii as outlined in Sec. |3.2| respec- 
tively. The underlying calculations take into account the eccentric- 
ity of the stellar encounters using Eq. 



different initial disc radius distributions (see Sec. 3.2 1 at 
the present age of the Arches cluster. The two plots show 
very similar distributions of models I and O. As discussed 
at the beginning of this section the only significant differ- 
ence occurs at the highest masses (M > 40 M0). Their 
suppressed disc-mass loss for constant initial disc radii 
(black bars) in the two times sparser model O (Fig. [5}d) is 
attributed to the less efficient gravitational focusing. So 
we will focus on the best-fitting model I for now (Fig.lsk). 

If we assume that disc radii are positively correlated 
with stellar mass (grey bars) then stellar encounters 
would act entirely destructive on discs of the most mas- 
sive stars and least harmful at the lower mass end. How- 
ever, in our preferred scenario of uncorrelated disc radii 
and stellar masses, stars in the mass range 10 — 30 Mq 
would suffer least encounter-induced disc-mass loss while 
the lowest and highest mass stars would lose two to three 
times more disc material, respectively. This result is 
anal ogous to simulations of the intermediate-mass ONC 
(see Pfalzner et al. 2006) yet for the much denser and 
more massive Arches cluster the average disc-mass loss 
is lowest for roughly three times more massive stars (see 
their Fig. 2). 

The reason for this minimum is a mass-dependent bi- 
variate encounter history: For low-mass stars the mass- 




1.0 1.5 

radial position [pc] 

Figure 6. Average normalized disc-mass loss as a function of 
distance to the cluster centre after 2.5 Myr of simulation time for 
a) model I and b) model O. The three lines mark three different 
mass groups of stars: 7 — 120 Mq (solid), 0.5 — 7Mq (dashed), 
and 0.1 — 0.5 Mq (dotted). The underlying calculations take into 
account the eccentricity of the stellar encounters using Eq. m. 
Constan t init ial disc radii have been assumed for this calculation 
(see Sec. [3^. 

loss occurs either through a few strong parabolic en- 
counter events with high-mass stars or frequent weak 
hyperbolic encounters with low- and intermediate-mass 
stars. High-mass stars lose their discs via a steady nib- 
bling by many gravi tationally focused e ncounters with 
stars of lower mass (lOlczak et al. 2010). For stars in 
the intermediate mass regime both interaction channels 
are reduced such that encounter-induced disc-mass loss 
is least efficient. 

The distribution of the resultant disc-mass loss with 
distance from the centre of the Arches cluster at its 
present age is presented in Fig. l6k for model I. It con- 
firms the expected decrease of the average normalized 
disc-mass loss with distance from the cluster centre as a 
consequence of the decreasing local density. However, be- 
yond roughly two half-mass radii - the here used formal 
cluster boundary - the average normalized disc-mass loss 
rises again. This is the footprint of strong stellar inter- 
actions that occur preferentially in the central region of 
the cluster. These interactions increase simultaneously 
the normalized disc-mass loss and probability for a rapid 
ejection of the encountered star. The escape speed trans- 
lates into distance from the cluster and explains its cor- 
relation with the normalized disc-mass loss outside the 
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cluster boundary. This can be seen best for the group of 
high- mass stars (sohd hne): the peaks in the histogram 
at large distances trace single escapers with huge mass 
losses. 

Unlike in model I, the radial distribution of the aver- 
age normalized disc-mass loss in model O (Fig. [6b) does 
not show the characteristic increase outside of the clus- 
ter. This is because the sample is dominated by stars in 
the tidal tails of the cluster that do not require strong in- 
teractions to become unbound from the cluster and thus 
mostly do not suffer increased disc-mass loss before es- 
cape. Note that the lower disc-mass loss of the high-mass 
sample compared to model I is again a consequence of the 
uncompensated gravitational focusing in model O. 

In the previous plots we have shown the evolution of 
the averaged disc-mass loss binned over different quanti- 
ties. Next we present a more detailed picture by showing 
the evolution of the spatial distribution of the disc-mass 
loss on a star-by-star basis at 1 Myr, 2.5 Myr, and 6 Myr 
in Figs. [7] to |9j The bottom panel shows the full ex- 
tension of the cluster, the two smaller top panels are 
zooms of 15 pc box half-width (left) and 3.5 pc box half- 
width (right), respectively. The grey scale represents the 
normalized disc mass rhd- In agreement with Figs. |3] 
and |6] the fraction of disc-poor stars (dark symbols) in- 
creases towards the cluster centre and with time. How- 
ever, the most prominent features are the growing tidal 
tails formed by stars escaping from the cluster in the tidal 
field imposed by the Galactic potential. By the end of 
the calculation the entire structure extends over an arc 
of nearly 150 pc length (Figj9|. 

However, we show in Fig. 
to the deep gravitational potential of the Galactic cen- 
tre only about 5000 stars or 7 % of the entire population 
are part of the tidal tails after 6 Myr of dynamical evo- 
lution. Here we define stars outside a distance of 3.5 pc 
to the cluster centre as escapers or equivalently as tidal 
tail members. The fraction of disc-less stars in the tidal 
tails remains roughly constant over time at about 5 % 
and amounts to 267 stars at 6 Myr. The escape rate is 
not constant but has periodical local maxima that coin- 
cide with the cluster's pericentre passages. It's shape is 
similar for star-disc systems and disc-less stars. 

For a deeper understanding of the dynamical history of 
the disc-less tidal tail stars in the following we investigate 
the temporal occurrence of three specific types of encoun- 
ters for each star. The events are referenced to the time of 
periastron passage of a perturber. At time tmax a star ex- 
periences the largest normalized disc-mass loss Amd due 
to an encounter. Note that for two events with equal rel- 
ative disc-mass loss Sm^ the first one causes a larger nor- 
malized disc-mass loss Awd, as is evident from Eq. ([s]). 
Hence tmax marks the event with largest absolute loss of 
the normalized disc mass m^, which is not necessarily 
coincident with the strongest perturbation of the disc. 
An event at tdstr marks encounter-induced disc destruc- 
tion. According to our criterion given in Section 3.2 this 
is the first encounter that reduces the normalized disc 
mass rhd below 0.1. The encountered star is considered 
to be disc-less subsequently. Finally, the time iiast marks 
the last encounter before escape from the cluster. Note 
that this does not necessarily imply the event induced an 
ejection from the cluster. It is possible for a star to be 
accelerated above escape speed in a previous encounter 



that despite its proximity 



and to experience only a weak last encounter on its way 
out of the cluster. However, the low probability for such 
an event implies that in most cases the encounter at tiast 
is energetic enough to unbind the star from the cluster. 

In the following we use the coincidence of these three 
dynamical events to split the entire set of 267 disc-less 
tidal tail stars at the end of the sim ulation, T* , into four 
disjoint subsets (marked in Fig. 11 by the symbols given 
in parentheses): 

T' = {xeT'\ tiast = idstr = tmax} (black filled) , 
Tb ={x eT' \ tiast = idstr > ^max} (black Open) , 
T' = {x eT' \ tiast > idstr} (grey filled) , 

— {x ^T* \ tlast < idstr} ■ 

Subset contains all the stars that lost their disc in the 
most perturbing encounter, being simultaneously ejected 
from the cluster. This implies highly energetic interac- 
tions. The difference to subset is that disc destruction 
for those stars occurred in a subsequent weaker event 
before escape from the cluster. However, here the last 
encounter was probably still quite energetic as simulta- 
neous disc destruction and ejection are less probable for a 
weak interaction. In the third set Tq we group stars with 
disc destruction before the last encounter. The dynami- 
cal history of these stars is potentially quite diverse. The 
disc might be destroyed in a rather strong encounter that 
does not immediately lead to ejection from the cluster. 
This is either because the encounter is not sufficiently 
energetic to unbind a star residing deep in the cluster 
potential or the potential escaper is back-scattered in a 
subsequent encounter. Alternatively, disc destruction is 
driven by a series of weak encounters. Eventually the 
disc-less star either evaporates out of the cluster or is 
ejected via an energetic encounter. Membership in the 
last subset 7^ is restricted to stars with disc destruction 
after escape from the cluster. The cardinality of the four 
sets is \Tl\ = 88, ITb'I = 79, ITc'l = 100, and ITj^l = 0. 

We expect the energetic encounters of stars in set 7^ 
to involve preferentially a massive perturber and rapid 
ejection from the cluster with high velocity. We expect a 
similar scenario for members of 7^ but a larger fraction 
of lower perturber masses due to the lower normalized 
disc-mass loss before ejection. In contrast, most mem- 
bers of set T' probably do not escape directly after disc 
destruction but remain bound to the cluster for some 
time. We expect disc destruction by perturbers of even 
lower mass without significant acceleration. An interest- 
ing consequence of the emptiness of 7^ is that disc de- 
struction does not occur in the tidal tails within the first 
6 Myr of cluster evolution. In other words all disc-less 
stars in the tidal tails have lost their disc before escape. 

In Fig. [TT] we present the dynamical properties of the 
disc-less tiHal tail stars in the three non-empty subsets. 
Panel a) shows that disc-less tidal tail stars have mostly 
low mass Mi (with a median of 0.23 Mq) and lose their 
disc at a small distance rdstr to the cluster centre (with 
a median of i/2i?hm ~ 0.4 pc). Only one very massive 
star with Mi > 10 M© becomes a disc-less tidal tail 
star as a member of T^. However, the distributions 
of Ml and rdstr do not differ significantly between the 
three subsets. In contrast, and as expected, panel b) 
clearly shows different distributions of perturber masses 
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Figure 7. Snapshots of the stellar population with encounter-induced disc evolution in the Arches cluster for the orbital model O after 
1 Myr of dynamical evolution. The bottom panel shows the full extension of the cluster, the two smaller top panels are zooms of 15 pc 
box half-width (left) and 3.5 pc box half-width (right), respectively. The grey scale represents the normalized disc mass rhj; the origin of 
the coordinate system marks the cluster centre-of-mass. The symbols represent three different stellar mass ranges: 7 — 120 Mq (triangle), 
0.5 — 7M (7i ( rectangle), and 0.1 — 0.5 Mq (circle). The underlying calculations take into acco unt the eccentricity of the stellar encounters 
using Eq.^4j. Constant initial disc radii have been assumed for this calculation (see Sec. |3.2| l. 
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Figure 10. Time evolution of the number of star-disc systems 
(grey line, left scale) and disc-less stars (black line, right scale) in 
the tidal tails of model O. The underlying calculations take into 
account the eccentricity of the stellar encounters using Eq. 
Constan t init ial disc radii have been assumed for this calculation 
(see Sec. [Oil. 

Af2,dstr that caused disc destruction and stellar veloci- 
ties Udstr after disc destruction. The lowest perturber 
masses for the three subsets TJ, 7^, and T' are 3.4 M©, 
0.72 Mq, and 0.13 Mq, respectively. The velocity dis- 
tributions differ strongly between set T* and the union 
set Tab = Ta ^"^b- While half of the stars in Tq have 
'^dstr < 16km/s, this is the case for only 3 of 167 stars 
in Tab- formally most of these latter stars are esca- 
pers with v > 2cr^, where cr^ ~ lOkm/s is the three- 
dimensional velocity dispersion of model O. We thus con- 
clude that nearly two thirds of the disc-less tidal tail 
stars lose their disc in a strong scattering event with a 
massive perturber. However, the distribution of the nor- 
malized disc-mass loss Am^ as a function of the time 
delay At between disc destruction and escape in panel c) 
shows that not all of these formal escapers are imme- 
diately expelled from the cluster. A fraction of 39/167 
stars from the union set T^b remains in the cluster for 
more than 1 Myr, even after very strong encounters with 
Arhd > 0.9. So even the strongest perturbations do not 
necessarily lead to the ejection of stars from the cluster. 
On the other hand does a small fraction of 17/100 stars 
from set T* escape within a few 0.1 Myr. In total af- 
ter 6 Myr 145 of 267 disc-less stars have escaped within 
0.5 Myr to become part of the tidal tails. We find no 
correlation between the normalized disc-mass loss caus- 
ing disc destruction and the subsequent time to escape 
from the cluster. 

In agreement with the rapid ejection of more than half 
of the disc-less stars, Figs. [T] and [8] show that there is 
a clear difference in the spatial distribution of star-disc 
systems and disc-less stars in the tidal tails. The disc- 
poor stars are preferentially located at large distances 
perpendicular to the tidal arms (i.e. along the orbital 
motion), and close to the cluster along the direction of 
the tidal arms (i.e. radially to the Galactic Centre). This 
is a consequence of strong interactions that remove disc 
material and accelerate stars sufficiently to overcome the 
deep cluster potential and escape in arbitrary directions 
- not preferentially via the Lagrange points. The spe- 
cific angular momentum L of these stars is less altered 
than that of stars escaping radially through the Lagrange 
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Figure 11. Distributions of characteristic dynamical quantities of 
disc-less tidal tail stars from the three sets (black filled squares), 
7^ (black open squares) , and 7^ (grey filled squares) . The panels 
show a) stellar mass Mi vs. relative distance to the cluster centre 
after disc destruction r^st^ in terms of cluster half-mass radii r^i^ , 
b) perturber mass Af2,datr that caused disc destruction vs. stellar 
velocity Vdstr after disc destruction, and c) normalized disc-mass 
loss Arhd vs. time delay At between disc destruction and escape. 



points. We demonstrate this effect in Fig. 12 by plotting 
the relative specific angular momentum of tidal tail stars 
AL ~ L — Lc, where Lq is the specific angular momen- 
tum of the cluster centre-of-mass. Clearly, the modulus 
|AL| is lower for disc-poor stars and thus explains their 
smaller radial offset by Eq. (fTl). 
As a consequence of the different dynamical history of 
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Figure 12. Relative specific angular momentum AL of stars in 
the tidal tails of model O with respect to the cluster centre-of- 
mass as a function of their distance from the cluster centre at 
2.5 Myr. The grey scale represents the normalized disc mass m^. 
The underlying calculations take into account the eccentricity of 
the stellar encounters using Eq. Consta nt in itial disc radii have 
been assumed for this calculation (see Sec, \3.2f. 



tidal tail stars we observe the forraation of two classes 
of tidal tails in the present-day cluster (Fig. [s]): a 
shorter "disc-poor" tidal arm pair populated by stars 
expelled from the cluster and suffering huge disc-mass 
loss (> 80 %) and the "classical" arms that are formed 
by stars evaporating though the Lagrange points. This 
is how - according to our dynamical model - the Arches 
cluster would appear at present. 

With time the growing number of stars that are tidally 
stripped from the cluster mix with the disc-poor escapers 
and the two types of arms merge into one large structure 
(Fig. |9|. However, the disc-poor stars still populate two 
distinct regions: the edge close to the Galactic centre in 
the trailing arm and the far edge in the leading arm. 

5. SUMMARY AND DISCUSSION 

We have investigated the effect of stellar encounters on 
the evolution of protoplanetry discs using a numerical 
model of the Arches cluster, one of the most extreme 
Galactic star clusters observed so far. 

Our numerical method is based on a hybrid approach 
where we combine simulations of pure star cluster dy- 
namics and a parameter study of isolated star-disc en- 
counters]^ This approach involves necessarily some as- 
sumptions and simplifications that potentially mildly 
overestimate the effect of the encounter-induced pertur- 
bations of circumstellar material. Nonetheless, our work 
represents the most detailed study of the effect of stellar 
dynamics on the evolution of protoplanetary discs in a 
starburst cluster so far. 

We find that stellar interactions in the Arches clus- 
ter are dominated by hyperbolic fly-bys, with more than 
70% having eccentricities £ > 10. This is much differ- 
ent in sparser clusters wh ere nearly paraboli c encounters 
are by far most frequent ( Olczak et al.|2010 1 . According 
to Eq. ^ the normalized disc-mass loss in hyperbolic 
interactions with e > 10 is reduced by a factor < 0.38 

* "star-disc encounters" refer to the interaction of a star-disc 
system and a disc-less star and serve as a numerical simplification 
of the disc-disc encounters that occur in reality. 



compared to the parabolic case. Consequently, encoun- 
ters drive the destruction of only about one third of the 
stellar discs in the cluster core until its present-day age 
of ~2.5Myr. The process continues with reduced effi- 
ciency such that after 6 Myr half of the core population 
becomes disc-less. The extended period of disc destruc- 
tion in a starburst cluster is again in stark contrast to 
sparser clusters in which it ceases after roughly 1 My r 
due to significant cluster expansion ( [Qlczak et al.|2010| . 
Note that by neglecting primordial binaries in our mod- 
els we underestimate the efficiency of disc destruction for 
three reasons: (i) the larger number of stars increases the 
probability for encounters, (ii) binaries increase the prob- 
ability for strong few-body interactions, and (iii) discs are 
truncated in tight binaries. 

The inclusion of a Galactic tidal field in our model of 
the Arches cluster induces the growth of an extended 
pair of tidal tails. The fraction of disc-less stars in the 
tidal tails is nearly constant over time at ~5 % and hence 
much lower than in the cluster at ages later than 1 Myr. 
Only about half of the disc-less tidal tail stars have es- 
caped slowly via the Lagrange points while the other half 
has been ejected at high speed from the cluster centre 
directly after disc destruction. These stars are the dom- 
inant population of a weaker and shorter second pair of 
"disc-poor" tidal arms at the present age of the cluster. 
With time the disc-poor tails mix with the dominant 
"classical" tails and form new features: a concentration 
of disc-less stars at the inner edge of the trailing arm and 
the outer edge of the leading arm. 

In fact, we expect that dynamical features arising from 
strong encounters - such as the disc-poor tidal arms - 
could be even more prominent th an sh own by our model. 
For one, as outlined in Section |4.1[ our model of the 
Arches cluster in the Galactic tidal field is slightly un- 
derdense such that close encounters are underrepresented 
in the simulation. Second, we model a population of 
single stars only and ignore the dynamically important 
effect of binaries. Assuming the binary fraction in star- 
burst clusters is potentially as high as in clu stered star 
forming regions of lower mass (^50%: e.g. Beck et al 
2003[ [Kohler et al.|[2006^ , the probability for stronglew 



body encounters would incr ease sig nificantly ( |Hills|1975 
|He ggie|[T975 |Anosova]|1986[ ), boosting the production o: 



high-speed escapers and disc-less stars 

However, one could argue that taking into account ad- 
ditional disc destruction mechanisms that are indepen- 
dent of cluster dynamics could mask the difference be- 
tween the classical and disc-poor tidal arms. We discuss 
here shortly the potential contribution fro m two other 
important processes. For one, grain growth (van Boekel 
et al. 2003; Shuping et al. 2006; Brauer ct al. 2008; Blurn 
PrWurm||2008nBiFHstier et al.||20'l0p would hinder de 
tection 01 protoialanetary discs m the infrared, the only 
wavelength range available to detect signatures of indi- 
vidual circumstellar discs in the distant Arches cluster. 
Whether grain growth in its disc population is expected 
to be as efficient as in the less massive star-forming re- 
gions in our solar neighborho od depends on th e environ- 
mental conditions. Recently, Okuzumi et al. ( 2011a|b l 
have shown that dust growth m discs can be suppressed 
by significant ionization, a probable scenario given the 
strong X-ray irradiati on of the n umerous massive stars 
in the Arches cluster tCanto et al..2000: Oskinova 2005.). 
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Second, if the radiation field of UV and X-ray photons 
emitted by massive stars is strong enough to heat the 
discs efficiently, photoevaporati on will drive additional 
serious disc-mass loss (Johnstone et al. 1998 Adams 



et al. 2004; Clarke 2007j |Gorti et al.||2009[ [ 



However, 

Stolte et al. (2010) haveTound a fraction ot a few percent 
of A5V to BOV stars (corresponding to a mass range of 
2 — 20 Mq) in the Arches cluster that show indisputable 
features of protoplanetary discs. Regarding its current 
age of 2.5 Myr and enormous radiation field fed by more 
than 100 OB stars it seems that photoevaporation was 
not acting entirely destructive over this period. 

So we expect that characteristic features of encounter- 
induced disc-mass loss - such as the disc-poor tidal arms 
- could be still apparent at the present age of the Arches 
cluster. As the velocities of tidal tail stars do not dif- 
fer by more than 2(7^ « 20km/s from the large space 
motion of the cluster of ^200km/s (see also Fig. |12[), 
kinematically they c an be well separat ed from the field 
star population (e.g. Stolte et al.||2008 ). This makes the 
Arches cluster a highly interesting observational target 
to study the combined and distinct effect of photoevapo- 
ration and encounters on the evolution of protoplanetary 
discs in an extreme environment. 

Note that observations of the stellar population and 
its circumstellar discs in this dense and distant cluster 
require instruments with very high angular resolution 
9 < 0.1". Unlike successfully performed for the nearby, 
less dense young star clusters in the solar neighbourhood 



(e.g. Fcigelson & Montmerle 1999 Kastner et al. 
Preibisch et al. 2005| |Albacete Colombo et al.|!i 



2005 
WangI 



et al.|2009| ), members of the young Arches cluster cannot 



be separated from the older field stars using the Chan- 
dra X-ray Observatory because of its limited resolution 
^ ^ 0.5". Thus at present the only viable strategy to 
discern cluster members from field stars is kinematic se- 
lection based on long-term proper motion studies. 

A natural extension of such a study would include 
the ^--^4 Myr old, more extended Quintuplet cluster, that 
has been considered to be the Arches cluster's "older 
borther" (Figer et al. 1999). Indeed, core collapse of 



our cluster model at 2.5 Myr induces significant expan- 
sion (see Fig. 2]), naturally explaining the two starburst 
clusters as different evolutionary stages emerging from 
the same initial conditions. Hence the Quintuplet clus- 
ter represents the ideal observational target to verify the 
future evolution of the Arches cluster as predicted by our 
simulation. The characteristic concentration of disc-less 
stars at the edges of the extended tidal arms provides a 
critical test. 

We have thus initiated a multi-epoch K-hand cam- 
paign and L-band photometry of the Arches and the 
Quintuplet cluster to obtain a kinematically selected 
membership sample required to observationally detect 
the characteristic disc-poor tidal features. 

However, we note that near-infrared surveys are not 
ideal to trace the predicted spatial gradients in disk frac- 
tions. This is because this wavelength regime is domi- 
nated by emission from the inner disc regions (at <1 AU 
from the central star) while the encounter-induced disc- 
mass loss studied here mostly removes material from the 
outer disc parts (at ^100 AU). Hence, ideally, observa- 
tions would be carried out at millimeter wavelengths. 
The instrument of choice for this purpose is ALMA as it 



provides the high spatial resolution {9 « 0.1") required 
to resolve individual sources in the starburst clusters at 
the Galactic Centre. 
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APPENDIX 
ENCOUNTER TRACKING 
Numerical method 

The essential task of the encounter tracking software is to determine the strongest perturber of each star at a given 
time step, i.e. the body with the maximum gravitational force. Numerically this algorithm is solved by an effective 
oct-tree search (Kennel|2004). This scheme, repeated for each temporal data snapshot, provides the interaction period 



and the masses, positions and velocities of the perturber and the perturbed star at their pericentre passage. While the 
interaction period can be easily obtained by recording the time step of the beginning and the end of the interaction 
with the same perturber, the determination of the orbital parameters at periastron is more challenging. 

However, for the calculation of the disc-mass loss of a star with mass mi due to the fly-by of a perturber with 
mass m2, the determination of their minimum distance - the separation at periastron r-p - is crucial. Due to the 
discreteness of time in numerical simulations the closest approach of two particles Tmin at a time step ^min is usually 



not the minimum distance of their orbits r-p as shown in Fig. 13 As long as the interaction is weak and orbital 
periods are long the difference between the numerical minimum separation and the theoretical distance at periastron 
is negligible. However, in case of strong interactions differences can be quite large - large enough that the estimated 
normalized disc-mass loss differs by more than a tenth or the orbit may be even classified as non-regular (with respect 
to the criteria discussed below) and rejected. The solution is to interpol ate t he periastron distance r^ at time tp from 
the interaction parameters at some other time step as described in Sec. |A.2[ 

Once the full interaction parameters are determined the orbit is checked for applicability of Eq. ([s]) to calculate the 
disc-mass loss. This involves two criteria: the orbit (i) must be nearly Keplerian, and (ii) must represent a significant 
part of the full Keplerian orbit. The first criterion is satisfied if the eccentricities at the beginning, pericentre, and 
the end of the interaction do not differ by more than 10 % and if the orbit has positive curvature (i.e. it is concave 
with respect to the perturbed star). The second criterion requires the opening angle of the orbit section to amount 
to at least 10 % of the maximum opening angle of the corresponding conic section. Orbits that fulfil these criteria, 
termed "regular" in the following, are stored in a list that forms the backbone for the subsequent determination of the 
disc-mass loss. 



Periastron interpolation 

The interpolation of the periastron rp requires the relative position and velocity vectors of a perturber, r and u, at 
a time step t. We use the time step of the closest approach tmin and assume that the orbit is (nearly) Keplerian, i.e. 
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Figure 13. Illustration of the geometrical parameters of the periastron interpolation scheme. The black circle marks the position of the 
perturbed star, the grey circles mark the positions of the perturber at sequential time steps. The lower part represents an enlarged view 
and contains all quantities that are used for the interpolation scheme (see text). 



there is no significant perturbation by a third body. This assumption is of course not always valid in the dense cluster 
regions, where encounters within small A'^-body systems can occur. However, since the derived disc-mass loss is based 
on a parameter study that was restricted to single perturbers, contributions from higher-order perturbations can not 
be quantified. So the disc-mass loss can be underestimated in this case. 

Defining the total mass M = mi + m2 and reduced mass = 17111712/ M, we use the six constants of motion given 
by the specific angular momentum h = fx v and the Laplace- Runge-Lenz vector 

. r GM^ 

A = v X h r, 

r 

where G is the gravitational constant, to estimate the periastron distance rp and eccentricity e of the orbit from rniin 
and Vmin, 

CM 

The interpolation of the time of periastron passage 

where At is the time period it takes to move from r,„i„ to rp over the angle = Z(rmin, ^p), is solved in polar 
coordinates {r, ip). Again, we make use of the Laplace-Runge-Lenz vector to determine A</j via 



Ar = \ A\\r\ cosip . 



From conservation of angular momentum, 
and the polar equation of a conic section, 



(1+£>P 

1 + £ COS ip 
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one finds via separation of variables: 



dt = 



At = 



dip 



I (l + e cos ifiY 



I 



2 fAip 



The result of the integration depends on the eccentricity s, 

e sin Ay 



At 



4(1 + ef 



(1 + ecosi^)' 



2 tan-1 



tan 



tan 



3 Aip 



6 sin A(^ 



In 



(e-l)tan ^+Vi^ 



(£2-i)(i+ecosA<^) e^-ly/^2—i tan 4^ -Ve^ 



ife < 1, 
if e = 1, 
if e > 1. 



Finally we need to determine whether fmin was passed before or after fp in time, i.e. the sign of At, which is given by 
the following relation: 

<0 A^G(0,7r) ^ At>0, 

> ^ A(/? e (tt, 27r) ^ At<0. 



inin«^min 



